Analytic description of elastic electron-atom scattering in an elliptically polarized 
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An analytic description of laser-assisted electron-atom scattering (LAES) in an elliptically polar- 
ized field is presented using time-dependent effective range (TDER) theory to treat both electron- 
laser and electron-atom interactions non-perturbatively. Closed-form formulas describing plateau 
features in LAES spectra are derived quantum mechanically in the low- frequency limit. These for- 
mulas provide an analytic explanation for key features of the LAES differential cross section. For 
the low-energy region of the LAES spectrum, our result generalizes the Kroll- Watson formula to 
the case of elliptic polarization. For the high-energy (rescattering) plateau in the LAES spectrum, 
our result generalizes prior results for a linearly polarized field valid for the high-energy end of the 
rescattering plateau [A. V. Flegel et at, J. Phys. B 42, 241002 (2009)] and confirms the factorization 
of the LAES cross section into three factors: two field-free elastic electron-atom scattering cross sec- 
tions (with laser-modified momenta) and a laser field-dependent factor (insensitive to the scattering 
potential) describing the laser-driven motion of the electron in the elliptically polarized field. We 
present also approximate analytic expressions for the exact TDER LAES amplitude that are valid 
over the entire rescattering plateau and reduce to the three-factor form in the plateau cutoff region. 
The theory is illustrated for the cases of e-H scattering in a C02-laser field and e-F scattering in a 
mid-infrared laser field of wavelength A = 3.5 /jm, for which the analytic results are shown to be in 
good agreement with exact numerical TDER results. 

PACS numbers: 34.80.Qb, 34.50.Rk, 03.65.Nk 



I. INTRODUCTION 



The interaction of an intense laser field with atoms 
or molecules results in highly nonlinear processes whose 
spectra are characterized by plateau-like structures, i.e. 
by a nearly constant dependence of the cross sections on 
the number n of absorbed photons over a wide interval 
of n. These plateaus are well known for spectra of above- 
threshold ionization (ATI) and high-order harmonic gen- 
eration (HHG) [Trill- The rescattering picture 0-Q pro- 
vides a transparent physical explanation for the appear- 
ance of plateau structures: an intense oscillating laser 
field returns ionized electrons back to the parent ion, 
whereupon they either gain additional energy from the 
laser field during laser-assisted collisional events, thereby 
forming the high-energy plateau in ATI spectra, or re- 
combine with the parent ion, emitting high-order har- 
monic photons. High-energy plateaus originating from 
laser-driven electron rescattering were predicted also for 
laser-assisted radiative electron-ion recombination or at- 
tachment Qj 3 an d laser-assisted electron-atom scatter- 
ing (LAES) [9|,[l0j|. F° r laser-induced bound-bound (as in 
HHG) and bound-free (ATI) transitions, rescattering ef- 
fects are suppressed for an elliptically polarized laser field 
and completely disappear for circular polarization. In 
contrast, for laser-assisted collisional processes (such as 
LAES) a rescattering plateau exists even for a circularly 
polarized laser field [ll[ (cf. also Ref. [l2| )■ The classical 
rescattering scenario used to explain plateaus in LAES 
spectra for a linearly polarized field has been justified 



by a quantum-mechanically derived analytic formula for 
the LAES differential cross section which provides 
the rescattering correction to the well-known Bunkin- 
Fedorov [l4[ and Kroll- Watson 15 1 results. This formula 
factorizes the LAES cross section into the product of two 
field-free cross sections for elastic electron-atom scatter- 
ing with laser-modified momenta and a "propagation" 
factor (insensitive to atomic parameters) describing the 
laser-driven motion of the electron along a closed clas- 
sical trajectory. These three factors provide closed-form 
quantum expressions for each of the three steps of clas- 
sical rescattering scenario for the LAES process. 

Besides its fundamental interest for understanding bet- 
ter the physics of nonlinear phenomena, factorization 
of the outcomes for nonlinear laser-atom processes in 
terms of laser-dependent factors and factors describing 
the field-free atomic dynamics provides an efficient means 
for retrieving these atomic factors from measured spec- 
tra of strong-field processes. At present, such factoriza- 
tions form the basis for HHG and ATI spectroscopies 
that allow the retrieval of the photoionization cross sec- 
tions for the outer electron shells of atoms or molecules 
(from HHG spectra) (cf, e.g., Ref. [l6j) and differen- 
tial cross sections of elastic electron scattering from the 
positive ion of a target (from ATI spectra) (cf., e.g., 
Refs. [UGH). The factorization of HHG and ATI yields 
was first postulated based on numerical solutions of the 
time-dependent Schrddinger equation [19( (cf. also the re- 
view [2(|) and was then justified theoretically [within the 
time-dependent effective range (TDER) theory (2j l22j ] 
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for the case of a monochromatic field in Refs. 0, [24| for 
HHG and in Ref. [25| fo r ATI, and for the case of a short 
laser pulse in Refs. [H (for HHG) and 27] (for ATI). 
We note that in all the aforementioned studies only lin- 
early polarized laser fields were considered, in which case 
the theoretical treatment is simplified (due to the one- 
dimensional laser-driven propagation of the active elec- 
tron along the direction of laser polarization). However, 
although the driving laser ellipticity provides an addi- 
tional control parameter for intense laser-atom interac- 
tions, at present there does not exist a convincing justifi- 
cation for the factorization of the rates or cross sections 
of nonlinear phenomena in an elliptically polarized field, 
neither for laser-induced nor for laser-assisted processes. 

In this paper we show analytically that the LAES cross 
section in the region of the rescattering plateau cutoff 
may be expressed in factorized form (as the product of 
three factors) for the general case of an elliptically po- 
larized laser field. This result generalizes that for the 
case of linear polarization [l3[ and presents a rare ex- 
ample of a strong field process whose yield may be fac- 
torized for the case of a nonzero driving laser elliptic- 
ity. The results presented are obtained taking into ac- 
count the rescattering effects non-perturbatively within 
the TDER theory for collision problems [28] as reformu- 
lated for the case of LAES in a low-frequency, elliptically 
polarized field. Based on a detailed analysis of the two- 
dimensional closed classical trajectories of an electron in 
the laser polarization plane, we have obtained also an an- 
alytic estimate for the (non-factorized) LAES amplitude 
that describes the entire energy region of the rescattering 
plateau. Our analytic results are in good agreement with 
exact numerical TDER results. 

The paper is organized as follows. In Sec. II we provide 
the basic results of the TDER theory for the scattering 
state of an electron as well as for the LAES amplitude 
in an elliptically polarized laser field. In Sec. Ill we de- 
velop a low-frequency expansion for the key ingredient of 
TDER theory: the periodic function of time, f P (t), that 
enters the TDER result for the scattering state. This 
expansion allows one to approximate the scattering state 
as a sum of two terms: a zero-order ( "Kroll- Watson" ) 
term and a rescattering correction, which is responsible 
for the high-energy plateau in the LAES spectrum. The 
low-energy part of the LAES spectrum, described by the 
Kroll- Watson term in the LAES amplitude, is consid- 
ered in Sec. IV, while in Sec. V we provide a detailed 
analysis of the LAES amplitude in the rescattering ap- 
proximation, i.e., including the rescattering correction. 
In Sec. VI we present the factorized (three-factor) form 
for the LAES cross section in the rescattering approxi- 
mation, compare the LAES spectra in this approximation 
with exact TDER results, and discuss the influence of the 
laser ellipticity on key features of LAES spectra. Some 
conclusions and perspectives for further use of the TDER 
theory for description of LAES in an elliptically polarized 
field are discussed briefly in Sec. VII. Finally, in two Ap- 
pendices we present an alternative representation for the 



TDER LAES amplitude that we use for the exact nu- 
merical calculations within the TDER theory (Appendix 
A) and a brief description of the uniform asymptotic ap- 
proximation of an integral involving a highly-oscillatory 
function (Appendix B). 



II. BASIC EQUATIONS OF THE TDER 
THEORY FOR LAES 

A. Formulation of the problem 

We consider the scattering of an incoming electron hav- 
ing momentum p and kinetic energy E = p 2 /(2m) on a 
target atom in the presence of a long laser pulse approx- 
imated by a monochromatic, elliptically polarized plane 
wave having intensity I and frequency lo. We assume 
that both the electron energy E and the laser photon 
energy hu are small compared to atomic excitation ener- 
gies and that the laser parameters / and to are such that 
laser excitation/ionization of atomic electrons is negligi- 
ble. Under these assumptions, the electron-atom interac- 
tion can be approximated by a short-range potential U(r) 
(that vanishes for r > r c ). Thus, the LAES process can 
be described as potential (elastic) electron scattering ac- 
companied by absorption or emission of n laser photons 
(with n m i n = —[E/(Huj)], where [x] is the integer part 
of x). Thus, the momentum (or energy) spectra of the 
scattered electrons (the LAES spectra) are characterized 
by momenta p„ and energies E n — p 2 l /(2m) — E + nHui. 

For the electron-laser interaction, we use the dipole 
approximation in the length gauge, 

V(r,t) = -evF(t), (1) 

were F(t) is the electric vector of the laser field, 

F(t) = FRe(ee- 4 "*) , e ■ e* = 1. (2) 

The complex polarization unit vector e in Eq. ^ is pa- 
rameterized as 

e + in\k x el 
e= 7=4== J , -1 < r? < 1, (3) 

V 1 + V 2 

where e is a unit vector along the major axis of the 
polarization ellipse, the unit vector k defines the laser 
propagation direction, and rj is the ellipticity. With the 
definition ([3]), the laser intensity does not depend on 77: 
/ = cF 2 /(8ir). Along with 77, the degrees of linear (I) 
and circular (£) polarization are convenient parameters 
for describing an elliptically polarized field: 

1 - V 2 r *, 277 

£ = e-e=- -, £ = wc-exe* =- ^. (4) 

1 + rf L 1 1 + rf w 

Note that the scalar product of the polarization vector e 
with a unit vector u, defined by the two spherical angles, 
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6 U and cf) n , as u = (ecos0 u + [k X e] sin </> u ) sin 9 U 
fccos8 u , is complex and can be parametrized as 



where 



E 



is the quasienergy and u p 



9v 



arg(u ■ ej 



u • e = |u • e|e^ u 
|u ■ e| — sin 



(5) 



tan f u — V tan i 



For an analytic non-perturbative account of both the 
electron-laser and the electron-atom interactions in elec- 
tron scattering assisted by a low-frequency elliptically 
polarized laser field, we adapt the TDER theory [28[ 
for LAES to the case of a low-frequency field. The 
atomic potential U(r) is assumed to support a single 
(negative ion) weakly-bound state ipKimi (r) with energy 
Eq = -fi 2 K 2 /(2m) (kt c 1) and angular momentum 
I. In particular, I = corresponds to electron scattering 
from hydrogen or an alkali atom, and I = 1 corresponds 
to a halogen atom target. 

The key idea of the TDER theory is the same as in 
effective range theory for two stationary potentials, U(r) 
and V(r), which exert their influence on the electron pre- 
dominantly in two essentially non-overlapping coordinate 
ranges [29]: U(r) is important primarily for r < r c , while 
a long-range, external-field potential V(r) is important 
primarily for r k . Thus, in the region r c < r <C ft -1 , 
the low-energy electron may be considered as virtually 
free. In this case, as in effective range theory for low- 
energy electron scattering [3(J, only a single parameter, 
the l-wave scattering phase 6i for the potential U(r), 
determines the Z-wave component of the exact scatter- 
ing state ipp(v) in the region r c < r <C min(K ,As ) 
(k = V2mE/h = p/h): 



Yj* (r)V P (r)dfi r ~ r- 1 - 1 + ■■■ +B l (E)(r l + •••), (6) 



where the factor Bi(E) involves the phase shift 8i(k) and 
can be approximated by two fundamental parameters of 
the effective range theory: the scattering length (a;) and 
the effective range {r{): 



(2l-l)\\{2l + l)\\Bi(E) 
- l +nk 2 /2. 



„2i + l 



cot 8i{k) 



-n 



(7) 

The boundary condition for ^p(r) at small r is the 
key equation that allows one to represent the scattering 
state Vp( r ) outside the potential U(r) (i.e., for r > r c ) 
in terms of the two parameters of the problem, a; and r; , 
independent of the shape of U(r). 



B. Scattering state of an electron in TDER theory 

We seek the laser-dressed scattering state, x P p (r, t), 
of an electron in the LAES process using the Flo- 
quet or quasienergy state (QES) representation (cf., e.g., 
Ref. [HI): 

p (r, t) = e-^/^r, t), $ p (r, t) = $ p (r, t + 2ir/u), 

(8) 



e 2 F 2 /(Amuj 2 ) is the ponderomotive (or quiver) energy. 
The QES wave function $ p (r, f) is a periodic solution of 
the time-dependent Schrodinger equation: 



lh §i + £+ £ A " U{r) ~ V{V ' t] ) <1>p( J '' ' 1 = "' ( !M 



Owing to the time dependence of $ p (r, t), the boundary 
condition for the Z-wave component of $ p (r, t) at small 
r > r c must be modified compared to Eq. (|6|) by introduc- 
ing some time-periodic functions (as was done similarly 
in TDER theory for bound states in an elliptically po- 
larized field (21, 22]). Since V(r,t) lacks axial symmetry 
in the case of an elliptically polarized field, the /-wave 
component of < I > p (r, t) depends in general on the angular 
momentum projection m;. However, for small r > r c the 
potentials U(r) and V(r, t) can be neglected in Eq. ([9]), 
so that the /-wave component of any time-periodic solu- 
tion of Eq. (|!|]) is independent of mi and may be written 



^(r^pM^O- 

^[akji(xkr) + b k yi{,x k r)\e~ 



ikujt 



(10) 



where xj, = y / 2m(e + khu)/h, ji and yi are the reg- 
ular and irregular spherical Bessel functions (behaving 



respectively as 



and 



0), and Ofc 



and 6fc are constants. Replacing ji(xk r ) an d yi(x k r ) in 
Eq. (|10p by their expansions for x^r <C 1, defining the 
factor Bi(e + khaj) as proportional to the coefficient ratio 

cik/bk, and introducing coefficients f^ mi \ in which the 
index mi labels the angular momentum projection onto 
Yi mi on the left of Eq. (fT0|) , we obtain a generalization of 
the boundary condition © for a time-dependent inter- 
action V(r, t): 

k 

+ Bi(e + khw){r l + ■ ■ ■ )]/('"*' > e -ik<"t 

= [r- l - 1 + --- + Bi(e)(r l + ---)]f^(t) 

[(2Z + 1)!!] 2 h dt Jp 

where the effective range parametrization ([7]) for Bi(e + 
kfvjj) was substituted on the left of the equality in 
Eq. (jlll) in order to obtain the final result summed over 
k on the right in terms of the time-periodic function 



(Imi) —ikut 



(12) 



The desired solution of the exact equation ([9]) for the 
scattering states has the following general form: 



$ p (r ) t)= Xp (r,t) + $W(r ) i), 



(13) 
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where the "scattered wave" <&p ^ (r, t) is an outgoing wave 
at r — > oo, while the "incident wave" Xp( r >i) is the QES 
wave function of a free electron with momentum p in the 
laser field (i.e., the time-periodic part of a Volkov wave 
function) , 



x P M) 



= e i[r.P(t)-S p (i)]/fi 



(14) 



where 



s P (t) 



[P 2 (r)/(2m)-e]dT 



eF(t) 



; 2 A 2 (r) 
2mc 2 



dr, (15) 



and P(£) = p — (e/c)A(i) is the electron's kinetic mo- 
mentum in the laser field F(t) with vector potential A(t), 
where F(i) = -c^dA^/dt. 

According to the TDER theory [28}, the function 

$p ^ (r, i) in the outer region, r > r c [in which the po- 
tential U(r) vanishes], can be expressed in terms of the 
retarded Green's function G(r, t;r,t') of a free electron 

in the laser field F(t) and involves the function fp lmi \t) 
in the boundary condition (fTTj) . [Indeed, upon neglecting 
U(r), any solution of Eq. ^ can be represented as a wave 
packet composed of wave functions for a free electron in 
the field F(t).] For G(r,t;r,t') we use the well-known 
Feynman form: 



G(r,t;r',i') = -9(t-t')- 



-i 3/2 



27ri^(* - t') 



x exp[i5(r, t; r', t')/H], 



(16) 



r — > V r . Equations for the unknown functions / p (t) 
complete the construction of the scattering state 3> P (r, t) 
in TDER theory [cf. Eqs. {T3J and ([15])]. These equa- 
tions are obtained by matching the Z-wave components 
of $ p (r,i) [which are different for different values of mi, 
as noted above and as is clear from the explicit repre- 
sentation (fT5)) for $>p C \r,t)] at small r to the prescribed 
boundary condition (fTTj) . Due to the term x p (r, t) in 
Eq. (fT3|) . the resulting equations comprise a system of 
21 + 1 coupled inhomogeneous integro-differential equa- 
tions for the functions fp mi \t), with mi = —/,•■• , i. 
Because the derivation and analysis of these equations 
involve the same steps for both I > and I = (differing 
only in the complexity of the analytical transformations) , 
for greater clarity, in the rest of this paper we provide 
analytical derivations only for the case of I = ( "s-wave 
scattering"). (For an analytical treatment of a similar, 
though homogeneous, system of equations in TDER the- 
ory for bound states with I > 0, see Refs. [2~H |22T|.) 



C. Exact TDER LAES amplitude and differential 
cross section for s-wave scattering 

If the potential U(r) supports only a single weakly- 
bound s-state so that only the phase shift <5o (A;) is non- 
zero, then Eqs. (fTTj) and (fTSf simplify as follows: 



$(-)(r,i) = -— f dt'e^-^^Ut') 



where 9(x) is Hcaviside function and S is the classical 
action for an electron in the laser field F(t): 



xG(r,W), 



(19) 



S{v,t;v',t') 



2(t-t') 



mu>* 



[F(i)-F(t')] 



2mc 2 



<I.. p (r./i- [l.+B (e)) f p (t) + i r -^f t fp(t), (20) 



A 2 (r)dr - - [r • A(t) - r' • A(t')] . (17) where / p (t) = / p 00) (t) and 



The behavior of <J> p (r, i) as r — > required by the 
condition (|lip [namely, the /-wave component of $ p (r,i) 
should involve a singular term ~ r~ l ~ 1 Yi m[ (f)] may be 
ensured by Z-fold differentiation of G(r, i; r, t') over r' fol- 
lowed by the substitution r' = 0. [From the explicit 
form (jTHJ) of G, such differentiation does not change the 
asymptotic behavior of $ p (r,i) for r — > oo.] As a result, 
in a way similar to that for the TDER treatment of a 
quasistationar y q uasienergy state with an initial angular 
momentum I [2l[ [22| , the general TDER expression for 
$ p c ' (r, t) can be written as follows [28[ : 



B (e) = — a a 1 + r me/h 2 . 



(21) 



To match the function $ p (r,i) [cf. Eqs. (US]), (JUJ)] to 
the r — > boundary condition (1201) . we extract from the 
integrand in Eq. (fl9|) a term proportional to the field- 
free Green's function Go(r, i; 0, t') [given by Eq. (fiTJI) with 
F(t) = 0]: 



* p '°>(r,i) = 



27T^ ; 



2 



dt'[e i < t - t '^ h fp(t') 



$p sc) (r,t) 



2Trh 2 



.i+i 



=-/ 



xG{r,t;0,t')-fp(t)G o (r,t;0,t')]+—fp(t). (22) 



x^(V r 0G(r,i;r',i')U o , 



(18) 



where the differential operator 3^,(V r ) is obtained from 
the solid harmonic 3^ t (r) [= r Y5^(r )] by the substitution 



The integral in Eq. ([221 is now regular at r = 0. Setting 
then r = in x p (r, t), comparing the result for <I> p (r,i) 
at small r with Eq. (|20p , and introducing the dimension- 
less time r = wi, we obtain an inhomogeneous integro- 
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differential equation for / p (r) = f p (t — t/uj): 

. rpmuj d 
B o(e)/ P (r) + i— T-/p(r) 



fi, dr * 
= «-'W+I[/ p (r)], 



(23) 



mijj f°° dx 



r e (i/K)[ex/w+S(T,r-x)] 



x f P (r 



x 



(24) 

S(r = 0,t = 



where S p (r) ee S p (t = r/w), S(r,r') 
r/w;r' =0,*' = r'/w). 

As is usual, the LAES amplitude A n (p, p„) is deter- 
mined by the asymptotic behavior of the wave function 

oo. For s-wave scattering, 

-i4>(r,t)/h 



<£ p sc> (r,t) in Eq. (|T5|) as r - 
this behavior has the form: 

1> p sc >(r,d 

Ikj-»1 



X(P,Pn 



where 



e r* 

0(r,t) = -r-A(t)+ / 
R(r,t) = r+-^F(i), 



e 2 ^ 2 (r) 

1 

2mc 2 



(25) 



dr. 



and the summation over n involves all open channels with 
exchange of n photons, for which E n = E+nhu > 0. The 
LAES amplitude A n (p,p n ) may be expressed in terms 
of / P (r), 



MP, Pn) = f e mr+lS - (T)/r 7 P (r)dr, 

and the differential LAES cross section is given by 
dcr n (p,p n ) Pn 



dft r 



|A»(p,Pn)| 



(26) 



(27) 



For F(i) = 0, the function / p (F(t) = 0; r) ee / (p) 
reduces to the amplitude >4(p) for field-free s-wave elastic 
electron scattering on the potential U(r) in the effective 
range approximation (in which k — p/h), 



fo( P ) = KA(p), A( P ) 



+ r k 2 /2 - ik 



(28) 



For F(t) ^ 0, the function / p (r) is a key object of TDER 
theory, since it contains complete information on the 
modification of the electron-atom interaction by an ellip- 
tically polarized laser field in all LAES channels. Numer- 
ical evaluation of / p (r) is done most easily by converting 
the integro-differential Eq. ([2"B"| to a set of inhomogeneous 
linear algebraic equations for the Fourier-coefficients 
of / p (r) [cf. Eq. ££!]) in Appendix E]. The LAES am- 
plitude An(p, p n ) is then expressed in terms of fk and 
generalized Bessel functions [cf. Eq. (|A6I) ]. As follows 
from the boundary condition (1111) that determines the 



QES $ p (r, t) at small r [where the potential U(r) is most 
important] , physically, the coefficients fk govern the pop- 
ulation of QES harmonics of the scattering state ^ p (r,t) 
with energies e + khw that arise as a result of atomic- 
potential-mediated exchange of k photons between the 
electron and the laser field at small r. 

For s-wave scattering, the numerical results in this pa- 
per, referred to as "exact TDER results," are obtained 
by numerical solution of Eq. (|A1[) . followed by evaluation 
of the amplitude A n (p,p n ) according to Eq. (|A6j) . For 
p-wave scattering (1 = 1), the Fourier coefficients /^ lm ' > 
(where m = 0, ±1) of the periodic function (|12[) satisfy 
the system of Eqs. ([A"8]) . (JX9J, while the LAES amplitude 
is given by Eq. (|A13jl . 

An analytic evaluation of the LAES amplitude 
A n {p,p n ) can be performed in the low-frequency limit, 
in which case the low-frequency expansion for the solu- 
tion / p (r) of Eq. (|2"3"|) can be obtained. 



III. LOW-FREQUENCY EXPANSION OF f p (r) 

Because the classical actions S p (t) and S(t,t') in the 
inhomogeneous and integral terms of Eq. (|23[) oscillate 
with large amplitudes (~ u p /u>) for the case of an intense 
low-frequency field F(f), we seek the solution / p (r) of 
Eq. (j2"3")) in the following form: 



/ p (r)= 5p (r) e -'r^'[ £ M-]/(M ; 



(29) 



where <? p (r) and £(t) are smooth functions satisfying 
respectively the requirements that |dg p /dr| <C u p /(huj) 
and that the upper bound of £ (r) is of the order of u p . 

Before proceeding with an iterative solution of 
Eq. (|23p . we analyze first the low- frequency limit of the 
integral term I[/ p (r)] defined in Eq. For u p » Hui, 

the dominant contribution to the integral (j2"4")l comes 
from the neighborhood of the singular point x = 0, while 
the contribution from the domain x > can be evaluated 
using the saddle point method. Thus we approximate the 
integral I [after substituting there Eq. (|29p] as a sum: 



Z[/ p (r)]«Z(% p (r)]+5>[/ p (r)], 



(30) 



where integrals 1^ and X s are evaluated respectively 
over the vicinity of x = and at the saddle points x = 
x s > 0. In order to evaluate the term X^°\ we neglect the 
action S(t, t — x) (which is of order x 3 when x — > 0) in 
the integrand of Eq. (|24[) and approximate the function 
/ p (r — x) for x< 1 as 

/p^-^/pWe"^ 1 )- 6 !^'. 
We thus obtain for I^°> the following result: 



^ (0) [/ P (r)]=/p(r) 



muj f°° dx 



ixE(r)j{hw) 



if p (T)^2m£(T)/h. 



(31) 
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The result for X s is obtained by substituting Eq. ([29]) into 
Eq. (|24l) , followed by evaluation of X by the saddle point 
method: 

e (i/h)[ex s /u+S(r,T— x„)] 

Ufpir)] = / p (r - x a ) — (32) 



where we have introduced the dimensionless function 

0{t,x), 



1 d 



d S(r, t') 



Z(t') 



T —T — X 



(33) 



and the quiver radius, «o = \ e \F/ (mw 2 ), for free-electron 
oscillations in the field F(t). The saddle points 
solutions of the equation: 



x)+£{t-x) = 0. 



(34) 



The results (J3TJ) and (32]) for the integral terms I (0) 
and X s allow us to develop an iterative procedure for 
the solution of Eq. ([23l) for / p (r) in the low-frequency 
limit. To do that, we note that the saddle point contri- 
butions, X s , to the integral X in the approximation (|30[) 
are proportional to the dimensional parameter a^ 1 = 
/(|e|F), while J(°) is proportional to ^/2m£(r)/fi, 



where £ (r) 



Thus, the ratio of terms X s to I^ - 1 is 



determined by a dimensionless factor ~ hjj/u p . There- 
fore, the iterative account of terms X s (which, as we will 
show below, describe the rescattering effects in LAES) is 
valid at the condition 



<C 1. 



(35) 



It is worthwhile to emphasize that, besides the frequency, 
the condition ([331) involves also the field amplitude F, so 
that the low- frequency expansion for the QES $ p (r,i) 
can be called also a "strong-field" expansion, since al- 
ready for fvjj/u p < 1 the perturbation theory (in laser- 
atom interaction) for <& p (r, t) becomes divergent 

A. The zero-order approximation for / p (t) 

To obtain the zero-order approximation in the param- 
eter ftw/u p for the function / p (r) [/ p (t) ~ fp\ r )}, we 
note that the strongly-oscillating exponential factor in 
Eq. (|29l) is determined by the inhomogeneous term of 
Eq. (|2"3")l [taking into account Eq. (T5])]. so that 



£(t) = £(°>(t) 



2m 



(36) 
(37) 



The pre-exponential factor gp^ (r) can be obtained from 
Eq. ([2"B"]) after substituting there / p (r) — !• f p \r), omit- 
ting then the differential term (~ udgp^/dr) and the 



saddle point contributions to the integral ([24]) [retaining 
only the first term in Eq. ([301), given by Eq. ([31])], The 
result for <? p ^ (r) is thus 



< } (r) 



B [fW(r)]-iP(T)/fi 



kA[P(t)], (38) 



where the second equality, obtained using Eq. ([2"T]) . gives 
the amplitude *4[P(t)] [cf. Eq. ([25]) ] for laser-free elastic 
electron-atom s-wave scattering in the effective range ap- 
proximation, as a function of the time-dependent kinetic 
momentum P(t). 



B. The first-order (rescattering) correction to 

fi 0) (r) 

The first-order iterative correction fp^ to the zero- 
order result /p°' satisfies the equation obtained by sub- 
stituting 

/p(r) = /( 0) (r) + /W(r) (39) 

into Eq. (23)): 

P [£ (1) (t)] - l^mfW^j/WM 
= ^X s [/(°)(r)], (40) 

S 

where / p X) is taken in the form (J29J) with £(r) = £ (1) (t) 
and 3 p (r) = ffp (t). In deriving Eq. (30]). the differ- 
ential term in Eq. ([23]) is evaluated as follows: we ne- 



glected the terms involving dg p ^ jdr and dg p ^ jdr and 
combined the result of taking the derivative of the expo- 
nential [see Eq. ([2"9"]) j with the term involving Bo(e) [see 
Eq. (SI)] to obtain B [£ (1) (r)]. Also, we used Eq. (31]) 
for ^"'[/^(t)]. The explicit form for I s [/^ 0) (r)] follows 
from Eqs. ([32]) and ([33]) taking into account Eqs. ([36]) 
and ([57]): 

X s [fi 0) (r)}=9 p 0) (T-x s )— — , (41) 



where 



tp(T,x) 



a [x 3 s /3(t, x s )] 1/2 
x[p — q(r, r — a;)] 2 



2mhw 



uj r e 

\ -F(t-x) ■ [q(T,T-x) -p] 



e 2 P 2 iw 
+ Q 2 (r-x,T)/a;}, 

Q(r,r')=q(r,r')--A(r), 



(42) 

(43) 
(44) 



q(r;r0 = £j;A(^ = eF(r)-F (r0- 



C T — T 



(J T — T 
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For the saddle point equation (|34|) we have the following 
explicit expression: 



P 2 (t-x s )-Q 2 (t-x s ,t) =0. 



(46) 



One sees from Eq. (|4"Tj) for the terms X s [fp] on the 
right-hand side of Eq. ([40]) that the oscillating (exponen- 
tial) terms of the function fp\r) are partly determined 
through the phase functions tp(r,x s ), which depend on 

the saddle points x s . Thus, the desired function /p (t) 
can be expressed as a sum, 



itp(T,x s )—iS p (T)/K 



(47) 



where we have introduced a set of functions g P}S corre- 
sponding to each saddle point x s . Substitution of the 
form (g7J) for f p 1] into Eq. flU gives the following equa- 
tion for the pre-exponential functions gp 

\h s {r)e lv(T ' Xs) =0, 



(i). 



fc,(r) = ( B Q [E^{r)] -pm£^{r) ) 5 £>(t) 
3 p 0) (t - x s ) 



(48) 



I 1/2 ' 



(49) 



where the set of functions £s replaces £^ in Eq. (|4Q|) . 
Comparison of the exponential factors in Eq. (|47j) with 
that in Eq. ([2"9"]) gives the following definition for ^^(r): 



fi 1 )(r)= e -fi W -^(r,x s )~5 p (r)A] 
= Q 2 (r,r-a; s )/(2m). 



(50) 



In order to proceed, we assume, that any two different 
solutions, x a and x s i , of Eq. (]4"6"]) do not merge with varia- 
tion of r and, moreover, they are such that the inequality, 



(It 



[<p{t,x„) - <p(t,x s >)] >u p /(tuv), 



(51) 



is fulfilled for the range of values of p and parameters 
of the field F(t) considered in this paper. [The validity 
of Eq. (f5Tj) can be justified by a numerical analysis of 
Eq. (|46|) (cf. Section IVAp .] Under this assumption, the 
exponential factors in Eq. (|48|) can be considered as quasi 
orthogonal functions in the following sense: 



ft s (T)e^ T ' x ')-^ T ' x «'>W < h s {r)dT , s^s' 



Therefore, without losing accuracy, we can consider only 
the trivial solution of Eq. (|48|) . h s (r) = 0, which from 
Eq. (|4T)]) gives a set of uncoupled equations for the func- 
tions gp, s {r). 

Finally, taking into account Eqs. (|38|) and (|50|) . the 
preexponential factors gp} s (r), that determine the first- 
order correction /p (t) in Eq. (|47p. can be expressed in 



terms of two field-free elastic scattering amplitudes ([28]) 
with different, field-dependent momenta: 

„(i) ( T ) = K A p i T - x s )]A[Q(t, t - x s )} ^ 



The most remarkable consequences of Eqs. ([3"5]) and 
([52]) are that (i) both results involve an exact (within ef- 
fective range theory), non-Born field- free scattering am- 
plitude A(p) with laser-modified momentum; and (ii) the 
result (I5"2"]) involves this amplitude twice. Fact (ii) allows 
us to call the approximate result ([47]) "the rescattering 
approximation." Thus the existence of laser-induced re- 
collisions in laser-assisted collision processes becomes ap- 
parent already on the level of the QES wave function 
$ p (r,i), in which the electron-atom dynamics and its 
modification by a strong laser field are completely de- 
scribed by the function f p {t). The low- frequency analysis 
of the exact TDER equation ([2"3"]) . presented in this Sec- 
tion, allows us to obtain analytic closed-form expressions 
for the LAES amplitude ([26]) corresponding to the zero- 
order [Eqs. ([37). ([38]) ] and rescattering [Eqs. ([iTjl. ([52]) ] 
approximations for f p (r) and, therefore, for the scatter- 
ing state $ p (r, t). 



IV. THE ZERO-ORDER (KROLL- WATSON) 
APPROXIMATION FOR THE LAES CROSS 
SECTION 

Using the zero-order approximation / p (r) w /p (t) 
where /p (t) is given by Eqs. ([3"T]) and (|38p ]. we obtain 



for the LAES amplitude ([26]) the expression: 

-4i 0) (P,Pn) = ^ r A[P{r)y nr+iA ^dr, (53) 



where 



A„(r) = [S p „(t) - Sp(r)]/H = pcos(r - ^ t ), 
le^ 

P 



' fc 2 |e-t|, <^ t =arg(e-t), t = p T 



P. 



and the scalar product (e-t) is defined in accordance with 
Eq. ([5J. For the more general case of Z-wave scattering, 
a low-frequency analysis of the TDER equations leads 
to the expression (|53[) for the scattering amplitude in 
which ^l[P(r)] is replaced by [P(t), P„(r)], where 
P„(r) = p„ - (e/c)A(r), 



-4 (0 (P„P/) 



(2/ + l)(fc i fc / )'7? i (cosg i/ ) 
-1/OJ + nk 2 /2-ik 2l+1 ' 



(54) 



= |Pi./l/^i 7 ? j( 3 ') * s a Legendre polynomial, and 
= Z(pj,p/). Later, we will omit the superscript 
(I) denoting the amplitude for field-free scattering (j54|) . 
A(pi,Pf) = (pi,Pf), bearing in mind that A(pi , p/ ) 
contains information about the spatial symmetry of the 
bound state supported by the scattering potential. Note 
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that the amplitude A^ (p) for elastic s-wave scattering 
in Eq (|28p is isotropic and depends only on the modulus 
of the initial momentum. Thus, if necessary, the differ- 
ence between ^4^ s ) (p) and A^ ^ (pi , p/ ) will be indicated 
by using a different number of arguments. 

It is important to note that the "instantaneous" am- 
plitude A[P(t), P„(r)] that replaces A[P(t)} in Eq. ||53J) 
is not an elastic scattering amplitude (since |P(t)| ^ 
|P„(r)|). For the case of linear polarization {1 = 1) 
Eq. corresponds to Eq. (5.16) in Ref. [Hj], which 
involves the T-matrix off the energy shell. For the case 
of elliptical polarization, results identical to Eq. (|53|) were 
obtained in Refs. [H, Sj- 

In the low-frequency limit (p> 1), the amplitude (|5"3")l 
can be evaluated analytically using uniform asymptotic 
approximation methods for integrals [H, (cf. Ap- 
pendix [5J : 



^(0) _ jn^ntpt 



A+J n {p)+A- 



\/p 2 - n ' 



(55) 



where J' n {x) is the derivative of the Bessel function J n (x), 
A ± = ^[A e i(r+)±A el {T-)], (56) 

and Aei{r±) = A[P(r±),P n (r±)], where r = r± are sad- 
dle points of the integrand in Eq. (|5"3")l that satisfy the 
equation 



/ cisin(r - ip t ) = n. 



(57) 



Because of the equality |P(r±)| = |P n (r±)|, A e i(r±) is 
the on-shell amplitude for elastic field-free scattering with 
laser-modified momenta. This modification serves to dis- 
place p and p n by the shift Ap± = (|e|/c)A(r±). For 
the classically allowed region \n\ < p, Eq. (|57) gives: 



7r n 
t± = ipt + — ± arccos — , 
2 p 



A P± 



m 



± £[k x t]y/ p 2 — n 2 



2|e-t| 2 

+n[2£(e • t)e + (1 - £)(t - (k ■ t)kj\ 



(58) 



(59) 



where the degrees of linear (£) and circular (£) polariza- 
tion are defined in Eq. (j4j. Note that for the case of 
critical geometry, when (e • t) — > (and thus p s» 0), 
the result ([55]) . based on a saddle point analysis of the 
integral (j53|) . is not applicable. 

The result (|53|) for and the corresponding cross 
section, 



Pn 
P 



A+J n {p)+A- 



ipJ'nip) 



V? 



(60) 



may be simplified and reduced to the well-known Kroll- 
Watson formula [l5[ for the following particular cases of 
the laser polarization and the scattering geometry: 



(i) For the case of linear polarization (I = 1), we have 
Ap + = Ap_ = Ap, where 



Ap = 

so that A+ = A e i(p 
while the cross section 
Watson result UM : 



-mhijjn 



(e-t) 



f Ap,p„ + Ap) and A- = 0, 
60)) reduces to the original Kroll- 



da. 



(KW) 



(P,Pn) 



Pn j2( ^g e ;(P,P ra ) 

J n\P jri ' 
p «"P„ 



(61) 



where da e i/dVL = \A e i\ is the exact cross section for 
field-free elastic scattering and P = p + Ap, P„ = 
p n + Ap. Note that the momentum shift Ap for the 
case of linear polarization remains real in the classically 
forbidden region |n| > p. 

(ii) For the cases of forward and backward scattering 
along the major axis of the polarization ellipse (p||p„||e), 
Ap± in Eq. (I5T)1) reduces as follows: 



Ap± = 



mhu> 



en ± 



1 



(62) 



The collinearity of the vectors p, p„, and t gives the fol- 
lowing relations: |P(r+)| - |P(r_)|, P(r±) = -P„(r T ). 
Thus, A- = 0, A+ = A e i{T+) = A e i{r-), and the LAES 
cross section is given by Eq. (|6"T|) with P = p + Ap± and 
P„ = p n + Ap± , where Ap± is given by Eq. (|6"2|) . (This 
result is the same using either Ap + or Ap_.) 

(iii) For forward or backward scattering in the polariza- 
tion plane for a circularly polarized field (£ = 0), Eq. ([59]) 
gives 



Ap± 



mhuj 



tn ± (\k x t]v. 



(63) 



With Ap± given by Eq. (|63|) . the same analysis as for 
case (ii) is then valid. 

Note that other analytic expressions for the scatter- 
ing amplitude (|53|) were obtained in Refs. [U The 
LAES amplitude in the low-frequency approximation in- 
troduced by Madsen and Taulbjerg [34| [labelled the 
"peaked impulse approximation" (PIA)] has a form sim- 
ilar to Eq. ([55]). but involves the Anger and Weber func- 
tions [cf. Eq. (|B7|) in Appendix [B]. In Fig. [T] we com- 
pare the PIA result of Ref. [34J with the analytic result 
(1551) . the integral expression (|53"|) (within the effective 
range approximation), and the exact TDER results. The 
effective range theory parameters are those for e-H scat- 
tering: \Eq\ = 0.755eV, k = 0.236a.u., a = 1.453k" 1 , 
and ro = 0.623k _1 . One sees in Fig. Q] that the zero- 
order approximation (|53[) for the LAES amplitude re- 
produces well the oscillation pattern in the LAES spec- 
trum. It follows from Eq. (|55|) that these oscillations are 
well approximated by the Bessel function and its deriva- 
tive; they originate from an interference of two classi- 
cal electron trajectories corresponding to two different 
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n 

FIG. 1. (Color online) Differential cross section 

da n (p, p n )/df2 Pn for laser-assisted s-wave e-H scattering in 
the polarization plane (p||e, p„ _L k) for a scattering angle 
9 = Z(p,p„) = 20° in a C0 2 -laser field with ftw = 0.117 eV 
(A = 10.6 /an) and intensity I = 2.5 x 10 11 W/cm 2 . The 
incident electron energy is E — 1.58 eV and n is the 
number of photons absorbed (n > 0) or emitted (n < 0). 
Results are shown for two laser polarizations: (a) Ellip- 
tical polarization, with 77 = +0.58 (I = 0.5); (b) circular 
polarization, with r\ — +1 (£ = 0). Circles: exact TDER 
results [cf. Eqs. (|A6|I . (|27[) ]; dashed lines: results using the 
approximate amplitude (|53[1 ; thick solid (red) lines: Eq. (|60[) ; 
thin solid (blue) lines: peaked impulse approximation (PI A) 
result of Ref . [341 ] . 



times of collision, r + and r_. In contrast, the result of 
Ref. [HI exhibits an additional sharp oscillatory structure 
for da n /dQ as a function of n that stems from proper- 
ties of the Weber function; they do not have any physical 
meaning. 

As may be seen from Eq. ((58]) . the two real saddle 
points t± coalesce at the cutoff of the classically allowed 
region (i.e., for n = p). In the classically forbidden re- 
gion (\n\ > p), the solutions of the saddle point equa- 
tion (|57[) and the corresponding momentum shifts (|59[) 
become complex, so we analytically continue the result 
(1551) to this case. However, the complex displacements of 
momenta in the elastic scattering amplitude may cause 
some non-physical features in the LAES cross section. 
Thus, for example, for electron scattering with absorp- 
tion or emission of \n\ > p » 1 laser photons, the con- 
dition P 2 (r±)/(2ra) = — \Eq\ may be satisfied for ap- 
propriate laser parameters and geometry of the incident 
and outgoing electrons. For such conditions, the ampli- 
tude A e i has a pole, which corresponds to some point 
t = r' p ' (or to more than one point) on the complex 
plain of r. The coalescence of one of the saddle points 
t± with the point r^ p ' leads to the appearance of a non- 
physical resonant-like enhancement of the LAES cross 
section. (This fact is exhibited most clearly for the case 
of forward scattering and circular polarization.) Thus, 
for the general case of elliptical polarization, the result 
(1551) has limited applicability in the classically forbid- 



FIG. 2. (Color online) The same as in Fig. [T] but for the 
cases of linear (£ = 1) and circular (£ — 0) polarization, and 
for a larger range of n > 0. Circles: exact TDER results 
[cf. Eqs. (|A6|) . (|27|) ]; solid lines: results using the analytic 
amplitude ((55j. 

den region. For this case, an alternative analytic result, 
suggested in Ref. [lij], is obtained within an additional 
weak-field approximation and, therefore, is not applica- 
ble for the description of strong laser field effects, such 
as the plateau structures in LAES spectra. 

In Fig. [5] we present LAES spectra for e-H scatter- 
ing in linearly and circularly polarized C02-laser fields. 
The field intensity, electron energy, and scattering ge- 
ometry are the same as in Fig. [T] For both of these 
two limiting cases of the laser polarization, I — 1 and 
I = 0, as well as for the general case of elliptical po- 
larization (0 ^ I ^ 1), the zero-order (Kroll- Watson) 
approximation (|55p for the LAES-amplitude does not de- 
scribe the high-energy part of the LAES spectra (i.e., the 
rescattering plateau) , for which a proper account of laser- 
induced electron re-scattering from the potential U (r) is 
required 0, [l(| ■ For the low-energy plateau, the result 
([55]) is in good agreement with the exact TDER results, 
except for the case of the critical geometry [for which 
e • (Pn — P) = 0], as exhibited, e.g., by the pronounced 
suppression of the LAES cross section within the Kroll- 
Watson approximation as compared to the exact result 
for n — 2 and £ — 1 (cf. Fig. [5]) . This discrepancy is due 
to the fact that the scattering angle 8 = 20° is close to 
the critical angle, 6 cr = 21.05°, for the channel n = 2. 

V. THE RESCATTERING APPROXIMATION 
FOR THE LAES AMPLITUDE 

The rescattering correction An^ to the zero-order re- 
sult ([55]) for the LAES amplitude, 

A n (p, Pn )^A^+A^, (64) 

follows upon substituting / p (r) = fp\r) [where fp\r) 
is given by Eqs. (gTJ) and ([52])] into Eq. $Zjo§ to obtain: 

4(1) - V.A (1) 
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A {1) 



1 



where we have defined t' s 
4> s {t) = <p(r, x s ) +nr - 



^ e i<t>s(T) dT 

A[P(r' s )]A[Q(r,r' s )} - = n ,(65) 



(a) 14 



t — x s and, using Eq. (fT5|) . 
(p„ - p) • F(r), (66) 



ihuj 2 



where the functions ip(r,x) and j3(r,x) are defined in 
Eqs. ([42]) and ([43]) respectively, and a; s = x s (t) is defined 
implicitly by Eq. (|46|) . 

For the case of /-wave scattering, our analysis of the 
rescattering correction to the LAES amplitude yields 
again an expression like (|65j) . but with the field- free s- 
wave scattering amplitude A(p) [cf. Eq. (|28)) ] in the inte- 
grand of (|d3|) replaced by -4(pi,p/) [cf. Eq. ((ST 

^[P(^)]^^[P(^),Q(r;,r)], 
A[Q(t,t' s )] ^ A[Q(r,r' s ),P n (r)}, 



where Q(r,r') is defined by Eqs. flU} and (|45|) . 

The dominant contributions to the integral (|65p come 
from the vicinity of the saddle points r — Tfc, which sat- 
isfy the equation: 



2mhu 



d<j> s 
~dV 



= P 2 (r fc ) - Q 2 (r fc , r fe - a:,) = 0. (67) 



[In deriving Eq. (|67[) , use has been made of the rela- 
tions nfojj = (p 2 — p 2 )/(2m), dF/dr = (w/c)A(r), and 
Eq. (146)) .] The saddle point equations (|46l) and ([67)) com- 
prise a system of coupled equations having a transparent 
physical meaning. Upon colliding with an atom at the 
time moment t' s k = — x s , the electron changes its mo- 
mentum p to a field-dependent "intermediate" momen- 
tum q(rfc, t' s k ), which ensures the condition for return of 
the electron by the laser field back to the atom at the 
time moment Tk followed by a rescattering. The set of 
points x s determines the excursion times of the return- 
ing electron along different closed classical trajectories, 
while Eqs. (1461) and (|67l) represent the energy conserva- 
tion laws at the times of the first and second collisions. 
The argument Q of the field-free amplitude A in Eq. ([65)) 
is the instantaneous kinetic momentum of the electron in 
the laser field in the "intermediate" state with canonical 
momentum q [cf. Eq. (j4"4")) ]. 

A. Analysis of the saddle-point equations 



To evaluate the integral in Eq. ([6"5")) , it is instructive to 
analyze first the solutions of the system of coupled sad- 
dle point equations (|46| and ([67)) . Using dimensionless 
quantities, this system may be represented as follows: 



7 2 - u 2 
< v 2 - 



where 7 = up/(\e\F), 7, 

^T,T')=Cq(T,T')/(Mn 



2(7-iv)-±m(ee- iT ) = 0, (68) 
2(7„ - v) ■ Im(e e~ iT ) = 0, (69) 

wp n /(|e|F), and v = 



...} 




/ / 
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FIG. 3. (Color online) (a) The solutions of Eqs. ([68} (dotted 
lines), (|69p (dashed lines), and (|7U|) (solid lines) for different 
values of y n , indicated in the figure near the corresponding 
curve, 7 = 0.6, 7||7 n ||e, and polarization r\ = 0.5. The black 
arrows show the direction of movement of the coalescing solu- 
tions of the coupled equations system (|68|l , (|69[) with increas- 
ing 7„. The corresponding coalescence points [the solutions 
of the system of Eqs. (|68[) . (|70|) ] are indicated by the black 
dots labeled by the numerals s = 1,2,3,4. (b) The classical 
closed trajectories of the electron in the polarization plane 
of the field F(t). Thick solid (black) line: the coalesced (ex- 
tremal) trajectories corresponding to the solutions s = 1,2,3,4 
of Eqs. ([681), ©■ Thin solid (red) and thin dashed (blue) 
lines: the short and long coalescing trajectories correspond- 
ing to the solutions of Eqs. ([68)) . ([69)) . 



Despite the fact that Eqs. ([68)) and ([69)) are very simi- 
lar, their solutions in the plane of variables r and r' (or, 
alternatively, r and x, where x — r — t') differ because 
of the different ranges of the parameters 7 and 7„. In- 
deed, rescattering effects become important in the region 
of the LAES spectrum where "direct" scattering is clas- 
sically forbidden, i.e., beyond the region of validity of 
the Kroll- Watson result, where 7™ > \/2(l + t) — 7 (for 
TnlMI^) [llll38) < On the other hand, rescattering effects 
are most pronounced for low incident electron energy, i.e., 
E < 2u p or 7 < 1. 

The numerical solutions of Eq. ([6"81 for 7 = 0.6 (7||e) 
and Eq. ([69)) for different values of 7„ (7n||e) for the 
case of elliptical polarization with rj = 0.5 is shown in 
Fig.EJa). Fig.[3la) illustrates the fact that, for the range 
of parameters considered, Eq. ()69l) has at most two real 
solutions on the trajectory x = x s (t) of each solution 
of Eq. ([68]) . With increasing 7„, the points Tj_ tend 

toward each other and coalesce at r = t s for 7„ = 7^ . 
For example, the point 1 (ri = 1.453, x\ = 4.764) in 



Fig. [3] corresponds to 7„ 



(i) 



jn ' = 1.982, while the point 

(2) _ 



2 (r 2 = 1.523, x 2 = 7.368) corresponds to 7^ = 1.682 

The coalescence of two real solutions, t_^' — = r s , 

at 7„ = and their disappearance for 7„ > means 
that the first derivative of </> s (t) has a local minimum 
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at r = t s , while r and a; vary along the trajectory of 
the solution x — x s (t). Thus the point r = t s , x — 
x s = x s (t s ) satisfies two coupled equations: Eq. (|68| and 
d 2 4> s /d,T 2 — 0. The latter equation may be written as: 

(u - 7n ) ■ Re (e e~" ) (t-t') + Q 2 - Q- Q' -j- = 0, (70) 

where the following notations have been used: 

Q = iy(T,r')+Im(ee- lT ), 

O! = v(r, r') + Im(ee-' 1T '), 
and where dr' /dr is determined implicitly by Eq. (J68J): 

rfr^ _ Q ■ Q! 

dr ~ (7 - v) ■ Re(e e~ lT ') (r - r') + Q' 2 ' 

As one sees in Fig.[3ta), the solution (r s , x s ) of the system 
of Eqs. (|68|) and (|70|) depends only weakly on 7„. 

The solutions (t s ,x s ) may be grouped in pairs, labeled 
by two consecutive (odd and even) integer subscripts s 
[with the solutions (r s ,x s ) enumerated in order of in- 
creasing values of x s , starting with s = 1] . Analysis of the 
system of Eqs. (fTDj) shows that the odd- and even- 

numbered solutions of each pair correspond respectively 

(s) 

to greater and smaller values of 7„ . Moreover, the first 
pair of solutions (i.e., s = 1,2) provide two limiting val- 
ues for 7^: for 7„ > = j n ,max the system ([55]), ([Ml) 
does not have real solutions [the derivative dcf> s (r) / dr as 
a function of r and s has a global minimum at the point 

(jx, Xi)], while the two saddle points t± do not coalesce 
(2) 

for 7„ < 7„ . All other solutions (r s ,x s ) correspond 
to intermediate values of 7« . A similar alternation of 

(s) 

7„ with increasing x s exists also in the analysis of the 
ATI process and was described within the semiclassical 
rescattering model in Ref. [2[. 

Considering the classical motion of the electron in the 
laser field F(i) described by Newton's equation, mr = 
— eF(i), a closed classical trajectory may be found for 
each solution of the saddle point equations (|68p and (1551) . 
For the geometry 7n||7||e and an elliptically polarized 
laser field, these trajectories lie in the polarization plane 
(r = rnc + rj_[/c X e]) and are shown in Fig.^b) for dif- 
ferent values 7„. The two different rescattering times, 
t! and correspond to the long and short trajecto- 
ries respectively, while the coalescence point (t s ,x s ) cor- 

(s) 

responds to the extreme trajectory with j n = 7„ . The 
smallest value of x s (i.e. , x% ) is the return time of the elec- 
tron along the shortest extreme closed path. During its 
motion along this shortest trajectory, the electron gains 
the maximal classical kinetic energy -E„ im a X = 2u P 7^ max . 

With increasing s (for x s 1), the solutions r s tend 
to a constant value (independent of s), while the sets of 
solutions x s with odd and even s become equidistant: 
{x s+ 2 — Xs) — > 2tt. This fact is easily verified by con- 
sidering the solution of Eqs. ([55)1 and (J7UJ) in the limit 
x = t — t' ^> 1. For this case, assuming I7 • e| ^ and 



|7 n • e| 7^ 0, the system (|68|) . (1701 reduces to the much 
simpler system, 

7 2 + 27 • Im(e e - l{T - x) ) = 0, 7„ • Re(e e~ lr ) = 0, 

which has the following solution: 

t = V^ + \, (71) 

X2k-(i±i)/2 = <P ln ~ + 27rfc ± arccos f , (72) 

2|7-e| 

where tp^ n — arg(7„ • e), ip~f — arg(7 ■ e). The approxi- 
mate results ([7T|) . (fT2"j) are in reasonable agreement with 
the numerical solutions of Eqs. (l68l) . (|70l) beginning from 
the third pair of points (t s ,x s ) (for the example presented 
in Fig. [3l the relative error for T3 and X3 is less than 3% 
and 1% respectively, while for T4 and 24 the error is less 
than 2% and 0.6%). 

Finally, we note that the solutions (t s ,x s ) with even 
s do not contribute to the high-energy region near the 
rescattering plateau cutoff, while they are important for 
the low-energy part of the rescattering plateau. The 
boundary energy, E ni between these two regions of the 
LAES spectrum is governed by the parameter 7 n) which 

is the limiting value of jn as s — > 00, where ^ k for 
odd s approaches 7 n from above, while for even s 

approaches it from below. The equation for 7 y n follows 
from Eq. (|M|) : 7^ = 2|7„ • e|. Using the parametrization 
([5]) for the scalar product (7^ ■ e), the boundary energy 
E n = Ziip 7 ) 2 can be expressed as follows: 

E n = 4M p sin6»p, l (l +£cos20 p „), (73) 

where 9 Pn and <p Pn are the polar and azimuthal angles 
for the vector p„ (or 7„) in the basis (e, [k x e),k). 



B. Analytic formulas for the scattering amplitude 

(s) 

Due to the coalescence of the two saddle points 
for each s, the ordinary saddle point method must be 
modified in order to evaluate analytically the integral in 
Eq. (|65[) (which determines the LAES amplitude within 
the rescattering approximation) . For this purpose we use 
the modification suggested in Ref. [39| and used recently 
to obtain factorized results for HHG [23( and ATI [25[ 
yields. This modification consists in approximating the 
phase factor <f> a (t) by a cubic polynomial in the neighbor- 
hood of the point r = r s , followed by removing from the 
integral (|55|) the slowly-oscillating pre-exponential factor 
at t — t s and extending the range of integration to ±00. 
The amplitude An can then be evaluated analytically 
in terms of an Airy function, Ai(x) [37]. The standard 
uniform approximation (in which one approximates the 
smooth pre-exponential factor by a linear function in the 
interval between the points r = t± ) [HI, H(|) gives ap- 
proximately the same accuracy of results, but leads to 
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cumbersome formulas, which are less suitable for further 
analyses and physical interpretations. 

As discussed above, the function s (r) is approximated 
as follows: 



MT) ~ ^ + 2^ (T - ^ 



+ 



3huj 



(r-r s ) 3 , 



(74) 



where t' s — t s — x s , <f> s = (j> s (r s ), and the dimensionless 
factor a s is proportional to the third derivative of <f> a (r) 
at r = t s , where in calculating this derivative one must 
take into account the r dependence of x s (r), defined im- 
plicitly by Eq. (jrJSj) . One obtains 



2(^-7n)-Im(ee- iTs )+Aa s , 



(75) 



where v s = v{t s ,t' s ) and 



Aa, 



dT 3 



(x s (t) - Xs)(y - u s y 



The explicit form of Aa s is cumbersome. It is not pre- 
sented here because numerical evaluation shows that it 
gives only a minor contribution to the final results. 

Evaluating now the integral (p5|) . we take into 
account that the amplitudes A[P(t' s ),Q(t' s ,t)] and 
A{Q(t, t' s ), P„ (t)] depend only weakly on r in the neigh- 
borhood of the saddle points r± [which satisfy Eqs. (|4l)|). 
(j6"7) ]. Thus the amplitude .A, evaluated at r = r s , can 
be replaced by the (on shell) amplitude A e i of field-free 
elastic electron scattering. The result for the LAES am- 
plitude An^ is: 



= — V^ s -4 ei (P (s) ,Q' (s) )A ; (Q (s) ,Pl s) ), (76) 

Ck 



where P( s > = P(t' s ), P [ n s) = P„(r s ), Q' (s) = Q(r s ',r s ), 
and Q< s) = Q(r s ,r^). The factors £> s in Eq. ^TBJ) are 
expressed in terms of the Airy function: 



D„ = \ — 

u„ 



1/3 i$ t 



Ai(C-) 



1/3 



(77) 



where /3 S = /3(r s ) is given by Eq. (|43p . and 



[(P.^) 2 -(Q(-)) 2 ]/(2m) 
Mp[a s (fej/wp) 2 ] 1 /3 



(78) 



The expression (|76p may be simplified after further 
analysis and some additional approximations. First, in 
accordance with the above analysis of the solutions of the 
saddle point equations, the sum over s in Eq. (f76|) can be 
split into separate sums over odd s and even s. The sum 
over even s contributes to the scattering amplitude only 
in the low-energy part of the rescattering plateau defined 
by E n < E n [cf. Eq. ([75])]. Second, the contribution of 
each succeeding term of the sum in Eq. (|76p decreases 



because the coefficient D s decreases as D s ~ xj 3 ^ 2 . Fur- 
thermore, each succeeding odd (s = 2k + 1) term con- 
tributes negligibly to the scattering amplitude in the re- 
gion 7„ > 7„ because the Airy function Ai(£ s ) decreases 
exponentially for £ s > —1.019. Thus we assume that the 
term with s = 1 gives the dominant contribution in the 
region of rescattering plateau cutoff, that the term with 
s = 2 contributes most to the region of the onset of the 
plateau, and that other terms (with higher s) give correc- 
tions in the intermediate region. Finally, the amplitude 
for field-free elastic scattering is a smooth function of 
its arguments and changes only slightly with respect to 
variations of s having the same parity, owing to the quasi- 
equidistant feature of the solutions (t s ,x s ) [cf. Eqs. (|7Tj) 
and ([72"jl ]. These considerations allow us to approximate 

the amplitude An^ by separating the summation over s 
in Eq. (|76p into two sums (over odd and even s) and by 
removing the slowly varying amplitudes A e i , evaluated at 
the proper momenta, from under each summation. Since 
the main contributions to the sum (|76|) are given by the 
first terms of the two separate summations (for odd and 
even s), we assume that the momenta are the correspond- 
ing instantaneous kinetic momenta, evaluated at the (di- 
mensionless) times (ti,t{) for the odd s sum: P = PW, 
P„ = P { n \ Q' = Q'W, Q = QW, and evaluated at the 
times (t2,T2) for the even s sum: P = P^ 2 ), P„ = Pn\ 
Q' = Q' (2) ,Q = Q (2) . The result is: 

A£Hp,p») = — [D^A e i(P, Q')Aei(Q, P„) 

+ D^Aei(P,Q')A el {Q,P n )], (79) 

where D<°> = EZo D ^+i, = EZi D 2k- 

The approximate result (|79|) [as well as the more ac- 
curate result l|76p] shows that the LAES amplitude with 
account of rescattering effects is given by a sum of factor- 
ized terms: all effects of the scattering potential U(r) are 
collected in the two exact amplitudes A e i for field-free 
elastic electron scattering, while the factors D s [defined 
by Eq. (|77p in terms of an Airy function] depend only on 
the laser parameters. Therefore, neither the scattering 
amplitude nor the LAES cross section can be factorized 
over the entire rescattering plateau region as a product 
of only two (laser and atomic) factors; however, such a 
factorization becomes possible in the high-energy part of 
the rescattering plateau, due to the negligible contribu- 
tion of the second term in Eq. (fT9"|) in this region. 



VI. FACTORIZATION OF THE LAES CROSS 
SECTION IN THE RESCATTERING PLATEAU 
REGION 

A. Three-step formula for the LAES cross section 

In the high-energy part of LAES spectrum, we can 
neglect the second term of Eq. (JTSJ) for the LAES ampli- 
tude in the rescattering approximation as well as the first 
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(Kroll- Watson) term in Eq. ([H3)l, Substituting Eq. ([79]) 
into Eq. ([27)1 , we obtain a factorized result for the LAES 
differential cross section in the high-energy region of the 
rescattering plateau: 



da ( n\p, Pn ) _ da el (P,Q') 



Pm 



dflQ> 



W(p,p n ) 



dg e ;(Q,P w ) 

drip 



where the factor W(p,p„), 



W(p,p„) = 



-4HE 



Z> 



2fe+l 



, (80) 



(81) 



depends on the momenta p and p„ of the incident and 
scattered electrons through the explicit dependence of 
the instantaneous momentum P„ [= p„ — eA(r^)/c] 
in the argument of the Airy function in Eq. (|77p . and 
through the implicit dependence of the times t s — 
r s (p,p„) and t' s — r^(p,p„) on the momenta p and p„. 
Since Eq. (|80[) was obtained as a simplified, low-frequency 
version of the exact quantum results for the scattering 
problem, its expression in terms of three factors provides 
a convincing quantum justification of the classical three- 
step rescattering scenario of the LAES process for the 
general case of an elliptically polarized laser field. 

The cross section da e i(P, Q')/dflQ' in Eq. (|80]) de- 
scribes the elastic scattering of an electron with initial 
momentum p from the potential U (r) at the time mo- 
ment t' = t[/uj. Since the collision takes place in the 
presence of a field F(t), this term involves (instead of the 
momentum p) the laser-modified instantaneous momen- 
tum P of the electron at the moment of collision. The 
scattering direction is given by the vector Q', which is de- 
termined only by the vector potential of the elliptically 
polarized laser field and has the sense of an intermediate 
"kinetic momentum" of the electron in an "intermediate" 
state, immediately after the elastic scattering event at the 
moment t'. From this state the electron starts to move in 
the laser field up to the moment of the second scattering 
(or rescattering) . The cross section da e i (P , Q')/ c££1q/ , in- 
volving the instantaneous momenta P and Q', describes 
elastic scattering (since |Q'| = |P|), while the initial mo- 
mentum p changes to q (|p| 7^ |q|). In order to ensure 
the condition for return of the electron back to the origin 
[where the magnitude of the potential U(r) is maximal] 
at the moment t, the vector q = q(ri,r{) depends on 
two times: the time t' of the first collision and the time 
t = n/u of rescattering. The result of the rescattering 
at the moment t is that the electron with the interme- 
diate momentum q rescatters along the direction of the 
final (detected) momentum p„. This event is described 
in Eq. (|5Uf by the cross section da e i(Q, P n )/dflp n for 
field-free elastic scattering with instantaneous momenta 
Q and P„ (where |P„| = |Q|). 

The key factor in the factorized cross section (|80| is 
the propagation factor W(p,p n ). This factor describes 
the motion of a free electron in the field F(t) for the 
time At = t — t' resulting in the change of its initial 



kinetic momentum P to P„. Indeed, as is seen from 
the explicit form for D s= 2k+i in Eq. (|77p. the expres- 
sion (j5Tj) for W(p,p„) does not involve any dependence 
on the potential U(r) and is determined completely by 
the free electron motion in the field F(t). Our numeri- 
cal analysis shows that the sum over k in Eq. (|81l) con- 
verges rapidly for arbitrary electron energy E n in the 
rescattering plateau region, so that only the first few 
terms in this sum over the saddle points contribute sig- 
nificantly. These terms effectively take into account both 
short and long closed trajectories of the electron in the 
laser field. These trajectories correspond to the two solu- 
tions, t± \ of the saddle point equations (1551 . whose 
interference causes the oscillatory features in the LAES 
spectra, which originate mathematically from the behav- 
ior of the Airy function Ai(£ s ). The times t s — t s /lo 
and At s = x s /uj, which govern the magnitude of D s 
in Eq. (|77p. are respectively the moment of rescattering 
and the excursion time for electron propagation along 
the closed trajectory corresponding to the extreme path 
for which the sth pair of short and long trajectories co- 
alesce [as shown in Fig. |3£b)]. The numerator of the 
Airy function argument £ s in Eq. ([75)) represents the dif- 
ference between the kinetic energy of the electron with 
final momentum p„ and the maximum classical energy, 
(Q( 5 )) 2 /(2m), that can be gained by an electron with ini- 
tial momentum p in the laser field before the rescattering 
event. 

The physical interpretation of Eq. ([50)1 is most clear if 
we limit ourselves to the case of the high-energy plateau 
cutoff region in the LAES spectrum, for which only the 
first term of the sum in Eq. (fSTj) dominates and the factor 
W involves only a single term, D\\ 



W(p,p„) 



Pn 

a\p 



(82) 



For the case of linear polarization (£ — 1), the factor- 
ization ([50)1 with W(p,p n ) given by Eq. (|%2"|) coincides 
with that obtained in Ref. [13[. The result (15^) takes 
into account only the return of the electron along the 
first pair of short and long closed classical trajectories in 
Fig.[3Jb), while the terms with k > 1 in the sum over k in 
Eq. (15Tj) determine the correction to the propagation fac- 
tor in Eq. (|52"j) due to electron returns along other "odd" 
(with s = 2k + 1) pairs of short and long trajectories 
[cf. Fig.EpD)]. 



B. Comparison with the exact TDER results 

In Figs. S] and [5] we compare exact TDER results for 
s-wave scattering (cf. Section 1 of Appendix [AJ with the 
low-frequency analytic results (for effective range theory 
parameters do and tq corresponding to the case of e-H 
scattering). One sees that the analytic result (|75|) for 
the scattering amplitude describes well the entire rescat- 
tering plateau region of the LAES spectra [we find that 
the simpler two-term result ([75)1 for An\p,p n ) provides 
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FIG. 4. (Color online) LAES differential cross section for 
forward s-wave e-H scattering (p||p„||e) as a function of the 
number n of absorbed laser photons in an elliptically polar- 
ized C02-laser field with huj — 0.117eV (A = 10.6 )im) and 
intensity / == 2.5 x 10 11 W/cm 2 for three different degrees of 
linear polarization I (= 0, 0.5, 1), and incident electron energy 
E = 1.58 eV. Thick solid lines: exact TDER results; dotted 
lines: result (|76p for the LAES amplitude; dashed lines: the 
three-step formula (|80[) : thin solid lines: Eq. (|80l) with ap- 
proximation (|82|) for the propagation factor. Vertical dotted- 
dashed lines mark the position of the boundary [cf. Eq. (|73|) ] 
between the two regions of the rescattering plateau. Arrows 
indicate the positions of the interference maxima and the 
plateau cutoffs. 



the same accuracy in describing the rescattering plateau] . 
For the high-energy part of the plateau (E n > E n ), the 
three-step formula (|80l) is in good agreement with the ex- 
act results. Moreover, the main contribution is given by 
the term corresponding to the shortest excursion time of 
the electron along the closed trajectory [cf. Eq. (jHH)]- The 
account of the longer trajectories [given by the terms in 
Eq. (|5TJ) with k > 0] provides a correction to the result 
([52]) in the spectral region between E n and the energy 
corresponding to the last (closest to the plateau cutoff) 
oscillatory minimum. 

Our analysis shows that the agreement between the an- 
alytic formula and the exact results in the cutoff region 
worsens for I — > 1 (cf. Fig. |4j. This fact is connected 
with the loss of the contributions to the scattering am- 
plitude of the intermediate QES-channels with negative 
quasienergies e n = E + u p + nhu [cf. Eq. (IA7|) ] when the 
saddle-point approximation for the exact TDER equa- 
tions was made. The effect of the closed channels on 
the LAES amplitude is not considered in this paper. We 
just note that the contributions of the closed channels to 
the LAES cross section in the high-energy plateau region 
noticeably depends on the laser intensity and the inci- 
dent electron energy E for a linearly polarized field and 
disappears for the case of the circular polarization. 

The comparison of our analytic results with exact 
TDER results for p-wave scattering (cf. Section 2 of Ap- 
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FIG. 5. (Color online) The same as in Fig. [4] but for scat- 
tering in the polarization plane (p||e, Pn-L«0 with £ = 0.5 
and three different scattering angles 6 = Z(p, p n ) (0°,±30°). 
Bottom panel: the LAES angle-integrated differential cross 
section over the "forward scattering" hemisphere (see text). 



pendix[3} is presented in Figs. [6] and [3 where the effec- 
tive range theory parameters are those for e-F scatter- 
ing: \E 1 = 3.401 eV, k = 0.500 a.u., a x = 0.827k" 3 , and 
ri = -4.417k. The intensity, / = 6.92 x 10 12 W/cm 2 , of a 
mid-infrared laser field with hoj = 0.354 eV (A = 3.5 /j,m) 
and the incident electron energy, E = 4.78 eV, are cho- 
sen so that the ratios u v j (fiuS) and E/(Huj) have the same 
values as for s-wave scattering in Figs. [4] and [5j One sees 
that the accuracy of the analytic result (|76[) for the scat- 
tering amplitude and of the three-step formula (f80| for 
the LAES cross section for p-wave scattering is as good 
as for the case of s-wave scattering (cf. Figs. [4] and [5]). 



C. Discussion 

The analytic results ([50)1 - (|82p allow one to explain all 
features of LAES spectra in the region of the rescattering 
plateau, as shown in Figs. 0] and [5] for s-wave (e-H) scat- 
tering and in Figs. [6] and [7] for p-wave (e-F) scattering. 
Moreover, in the case that the field-free cross sections 
da e i/d£l have a smooth energy dependence, these fea- 
tures are governed by the propagation factor W(p, p n ) 
and are insensitive to the details of the potential U(r). 
In particular, the position of the plateau cutoff, as well 
as the positions of the maxima and minima in the os- 
cillation pattern below the plateau cutoff, are described 
quantitatively with high accuracy by the properties of 
the Airy function Ai(£i) [where fi is defined in Eq. (|78j) ] 
(cf. similar analyses for high-energy HHG and ATI spec- 
tra in Refs. 23, 25]). If the energy difference in the nu- 
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FIG. 6. (Color online) The same as in Fig. [4] but for p-wave e- 
F scattering in a laser field with huj = 0.354 eV (A = 3.5 /-tm), 
I — 6.92 x 10 12 W/cm 2 , and incident electron energy E — 
4.78 eV. 
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FIG. 7. (Color online) The same as in Fig. [5] but for p-wave 
e-F scattering. The laser field parameters ui and /, and the 
electron energy E are the same as in Fig. [6] 



merator of £i in Eq. (1751) is positive, the Airy function 
(and hence the LAES cross section) decreases exponen- 
tially with increasing p n . In contrast, Ai(£i) oscillates 
for d < with the position of its first maximum at 
Ci = z\ = —1.019. This value of £i thus determines the 



position of the plateau cutoff (E n 



the LAES spectrum, where 7 C satisfies the transcendental 
equation obtained by equating £i to z\. 



v x + 2(7 C — Ui) ■ Im(ee 

1/3 ft , i n2/3 

z\ol{ (Tkjj/Up) , 



) 



- l = \ / 

/ ✓ 
/ ✓ 
/ ✓ 

-/ / 
Is 

e = o 


p — ft <; \ 
■>c-~ \ 

N \ 
\ \ 
\ \ 
\ \ 
\ > 
\ 

\ 

\ - 


= 30° / 

n _ f\o / \ ,.••** s' 
u — U X s 

^V0 = -3O° 







-45 -30 -15 15 30 0.2 0.4 0.6 0.8 1 

9 (deg) I 

FIG. 8. (Color online) The cutoff energy E c vs. scattering 
angle 9 for different values of the linear polarization degree 
I (left panel) and E c vs £ for different 9 (right panel). The 
scattering geometry is p||e, p n A-k and the laser parameters 
and the incident electron energy are the same as in Fig. f3] 



(83) 



where a± is given by Eq. (|75)l . v>i = v{t\,t[), and (n, t[) 
is the first solution (corresponding to the shortest re- 
turn time x\) of the system of equations (r70|) with 
In = 7c- In other words, the cutoff parameter 7 C is given 
by the joint solution of the coupled system of Eqs. (|68l) . 
([70]) . and (JS3J)- For an arbitrary ellipticity (including the 
case of circular polarization), an analytic expression for 
7 C may be found only as a polynomial interpolation of the 
exact numerical solution of Eqs. (|68|) . (f70]) . and (|83|) and, 
in general, this interpolation has a cumbersome form be- 
cause of its dependence on the many parameters of the 
problem (such as, e.g., the scattering geometry, the scat- 
tering angle, the ellipticity, the incident electron energy, 
and the laser intensity). Thus we show in Fig. [5] the nu- 
merical solutions of the transcendental equations for E c 
for scattering in the polarization plane for different values 
of the ellipticity and the scattering angle. 

As shown in Fig.[8j the cutoff position depends strongly 
on the scattering angle 9 for £ = 1: E c (9) w E c (0) — 
7.9u p 9 2 (cf. Ref. [H|]), while the angular dependence of 
E c (8) becomes smoother with decreasing linear polariza- 
tion degree £. For forward scattering along the direction 
of the major axis of the polarization ellipse, the depen- 
dence of E c (£) on £ is close to linear over a wide inter- 
val of incident electron energies E and laser intensities I 
[I = cF 2 /(8tt)]: E c (£)/u p « ai(E, F) + a 2 (E, F)£, where 
a\.2{E,F) are smooth functions of E and F (cf. Fig. [5]). 

Another noticeable effect seen in Fig. [5] is an asym- 
metry in the cutoff position with respect to the sign 
of the angle 9 for £ < 1 (cf. also Figs. [S] and [7J. (For 
the geometry p||e, p n _l_K, one has p n cos9 = p n • e and 
p n sin 9 = p n ■ [k x e] , so that the positive direction of 
9 coincides with the direction of the field rotation for 
r\ > 0.) This dichroic effect for the cutoff of the rescat- 
tering plateau in LAES spectra was predicted in Ref. [ll| ■ 

The oscillation pattern in the dependence of W(p, p n ) 
on p n originates from the interference of two classical 
electron trajectories, which merge at the cutoff with the 
shortest extremal trajectory and which were taken into 
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FIG. 9. (Color online) Scaled positions (E n ,msLx.) of the 
plateau cutoff and of the oscillatory maxima closest to the 
cutoff (marked by k = 1,2,3) vs. scaled laser intensity 
I/Io for forward scattering in a circularly polarized field 
and for two electron energies: E = 13.5ftw (solid lines) 
and E = 19.5ftu (dashed lines). [Scaled units of intensity, 
Io = 1.5 x 10 12 W/cm 2 , and of energy, \Eo\ = 0.755 eV, corre- 
spond to e-H scattering (see text). For the case of a CC>2-laser: 
E = 1.58 eV and 2.28 eV.] 

account in evaluating the LAES amplitude (cf. discussion 
in Sec. IV A[) , This interference explains the oscillatory 
patterns in the LAES spectra below the plateau cutoff 
(for d < Z\) 1 which are known from numerical calcula- 
tions (cf. Ref. [ill and Figs. [4|- FTP and were discussed in 
Refs. [13| and [40|]. [In Ref. [4(| the origin of the oscil- 
latory patterns as a consequence of the interference be- 
tween real electron trajectories was established by taking 
into account the scattering potential U(r) perturbatively 
within the strong-field and uniform approximations.] 

The positions of the minima/maxima of the interfer- 
ence oscillations may be found in the same way as for 
the cutoff position, i.e., by solving the system (TTTjl) . 
and Eq. (J53J) for j n = 7„, mi „/ max , replacing z\ by z k 
[where z k are the positions of the zeros and the maxima 
of Ai 2 (£i)]. For k > 2, the values of z k are well approx- 
imated by equating to irk/2 the argument of the sine 
function in the asymptotic form of Ai(— |£i|) for large 

idi na, 

Ai(-lCil) oc |Ci|- 1/4 dn (|lCi| 3/a + J) • 

The maxima/zeros of Ai 2 (£i) (and hence the max- 
ima/minima of da n /dQ) correspond to odd/even k in the 
relation 

C il = Zk = -0.25[27r(2fc - 1)] 2/3 , k > 2. 

The estimated positions of a few maxima in the LAES 
spectra closest to the cutoff are indicated in Figs.H]-[7]by 
arrows. One sees that these positions coincide well with 
the positions of the maxima in the exact TDER results. 

We have found that the positions of the maxima or 
minima in the oscillatory LAES spectra depend on the 
scattering angle and on the laser polarization in much 



the same way as shown for the cutoff position, E c (9,£), 
in Fig. [5] However, the distance between the positions 
of the maxima or minima for fixed £ and 9 depends es- 
sentially only on the laser intensity and scales as I 1 / 3 . 
This fact is shown in Fig. [9] for the case of forward e- 
H scattering in a circularly polarized (£ = 0) field for 
two values of the electron energy E. [The scaled unit 
of intensity, Iq, in Fig. [S] is denned as Iq = cFq/(8it), 
where F = y/8m\Eo^ / (\e\h) . Thus for e-H scattering 
(\E \ = 0.755 eV), I = 1.5 x 10 12 W/cm 2 .] Note that for 
a linearly polarized field, the same intensity dependence 
for the positions of the maxima and minima was found 
analytically for LAES [13] and for ATD [25[ processes. 

Because of the sensitivity of the oscillatory patterns 
in the LAES spectra to the scattering angle (cf. Figs. [5] 
and[7J), the angle-integrated spectra are smooth, as shown 
in the bottom panels in Figs. [3] and [71 in which the in- 
tegration was performed over the "forward scattering" 
hemisphere: < 9 Pn < 180°, -90° < <j) Pn < 90°, where 
0p n and (j) Pn are the polar and azimuthal angles for the 
vector p„. For this case, one sees in Figs. [5] and [7J that the 
simple analytic result ((80)) with propagation factor (|82j) 
provides good agreement with the exact TDER results 
over the entire rescattering plateau. 



VII. CONCLUSIONS AND PERSPECTIVES 

Nowadays the manifestation of field-free atomic dy- 
namics in strong field processes and the retrieval of in- 
formation on this dynamics from the measured outcomes 
of laser-atom interactions are attracting increasing inter- 
est. For HHG and ATI processes, this dynamical infor- 
mation can be obtained theoretically most convincingly 
using well-developed algorithms for direct numerical so- 
lution of the time-dependent Schrodinger equation. How- 
ever, for laser-assisted collisions, numerical algorithms 
for calculating the scattering state wave function in an 
intense, low-frequency laser field have not yet been de- 
veloped, even for the case of linear laser polarization. 
Moreover, the widely-used strong field approximation is 
not applicable for this purpose since for an electron in 
the continuum it treats the scattering potential pertur- 
batively, using the Born approximation. Thus for colli- 
sion problems, non-perturbative approximate theories or 
exactly-solvable models play an essential role in providing 
a deeper understanding of the influence of the scattering 
potential on laser-assisted collision processes. 

In this paper, we have obtained quantum-mechanically 
(in the low-frequency limit) analytic expressions for cross 
sections of electron scattering from a potential in the 
presence of an elliptically polarized laser field using 
TDER theory, which permits one to obtain not only an 
exact numerical solution for the LAES problem but also 
simple analytic results for a number of limiting cases. 
Our analytic derivations are based on the analytic repre- 
sentation of the exact TDER scattering state $ p (r, t) in 
Eq. (|13[) as a sum of two terms: the "zero-order" term, 
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which corresponds to the low-frequency, Kroll- Watson re- 
sult for the scattering state [cf. Eq. (5.12) in Ref. [l5j]. 
and the "rescattering correction," which takes into ac- 
count the strong laser field modifications of the electron 
interaction with the scattering potential U(r) beyond the 
Kroll- Watson approximation. Since the Kroll- Watson 
term in the LAES cross section decreases exponentially 
beyond the classically- allowed region (for high n), the 
rescattering correction becomes dominant there and de- 
scribes perfectly the rescattering plateau in the high- 
energy region of the LAES spectrum. The high accu- 
racy of our analytic approximations for the exact TDER 
LAES amplitude is demonstrated by comparison of ana- 
lytic and exact numerical TDER results for the ellipticity 
and angular dependences of LAES spectra for two differ- 
ent cases: s-wave scattering (corresponding to electron 
scattering from hydrogen or an alkali atom; cf. Figs. [4j [5] 
for e-H scattering) and p-wave scattering (corresponding 
to a halogen atom target; cf. Figs. [HO for e-F scattering). 

The key results of this paper are the expression ([76)) 
for the LAES amplitude in the rescattering approxima- 
tion and the three-step formula (l80|) for the LAES cross 
section. The factorized result (1501) describes well the 
high-energy part of the rescattering plateau, while the 
non-factorized LAES amplitude ([76")) [as well as the two- 
term result ([79])] describes the LAES spectrum over the 
entire rescattering plateau region (cf. Figs. @] - [TJ. Af- 
ter substituting Eq. ([52")) for the propagation factor, the 
formula (I80|) provides a generalization of the result for a 
linearly polarized laser field [l3| to the case of nonzero 
driving laser ellipticity. 

The major limitation of the TDER theory model is 
that it takes into account only a single partial- wave scat- 
tering phase (in a given /-wave channel) for the poten- 
tial U(r) [HI, SH, whereas the entire set of phase shifts 
should be taken into account in describing elastic electron 
scattering by a neutral atom. However, this deficiency is 
compensated by the very clear and physically transparent 
interpretation of our key results ([To]) and ([50)1 . Indeed, 
(i) the quantum-mechanically derived factorized formula 
([50)1 agrees completely with the semiclassical three-step 
rescattering scenario for the LAES process giving, in fact, 
a quantum "replica" (or quantum justification) of this 
scenario; (ii) the account of rescattering effects in our 
analysis was performed non-perturbatively in the poten- 
tial U(r), so that the results ([78")) and (150)) contain the 
exact (non-Born) amplitude and cross section for elastic 
electron scattering by the potential U(r) within the ef- 
fective range theory; and (iii) the factors D s [cf. Eq. (1771) ] 
in Eq. (1751) . as well as the propagation factor W(p,p„), 
do not involve any parameters of the potential U (r) and 
thus are valid for any atomic target. [In particular, our 
results for the s-wave and the p-wave scattering show 
that these factors do not depend on the spatial symme- 
try of a bound state (if it exists) in an atomic potential 
U(r).} Therefore, it is reasonable to expect that a gener- 
alization of Eqs. (J7BJ and ([501) beyond the TDER theory 
may be performed quite straightforwardly, i.e., replac- 



ing the field-free scattering amplitudes A e i in Eq. ([76)) 
and the TDER cross sections da e i /dQ, in Eq. ([50)) by the 
amplitudes and cross sections for elastic electron scatter- 
ing by a particular real atom obtained from either ex- 
perimental measurements or accurate theoretical calcu- 
lations. Similargeneralizations of factorized TDER re- 
sults for HHG [lj and ATI [IH] yields to the case of real 
atomic targets have been shown to provide fine agreement 
with results of accurate numerical solutions of the time- 
dependent Schrddinger equation for the plateau cutoff 
region in HHG and ATI spectra. For LAES, the afore- 
mentioned generalization allows one to extend the for- 
mulas (|76p and (J80j to the case of atomic targets (such 
as inert gases) which do not support a bound state of an 
attached electron (i.e., a negative ion) in spite of the fact 
that the description of LAES within the TDER theory 
presented in this paper is not applicable for such cases. 
The use of the results ([76")) and (|80p for such cases that 
go beyond the present TDER theory will be described in 
a separate publication. 

The results in this paper become inapplicable for reso- 
nant electron energies, E ss /iHuj— \Eq\ — u p , at which the 
electron may be temporarily captured in a bound state 
ipKlmi ( r ) 01 the potential U(r) by emitting \i photons [4l| . 
and for threshold energies, E = khus, k = 1,2,..., at 
which the LAES spectrum may be affected considerably 
by threshold phenomena, corresponding to the closing 
(or opening) of the channel for stimulated emission of k 
laser photons by the incident electron (28|. Since both 
resonant and threshold phenomena have a purely quan- 
tum origin, when the discreteness of the photon energy 
nhw is essential, these phenomena disappear in the low- 
frequency approximation (hoj — > 0) used in the present 
work. An analysis of resonant and threshold phenomena 
for the LAES process in an elliptically polarized laser 
field will be published elsewhere. 

Finally, we note that, even for the simplest geometry, 
p||e, the ellipticity r\ of the laser field affects significantly 
the angular distribution (AD) of scattered electrons as 
compared to the case of linear polarization, because it 
destroys the axial symmetry of the AD that exists for 
rj = with respect to the direction of e. In particular, the 
ADs for rj ^ differ substantially for 77 = ±|iy|, thus ex- 
hibiting an elliptic dichroism effect whose detailed study 
for both the low-energy and the rescattering regions of 
the LAES spectrum is now in progress. 
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Appendix A: The matrix form of the TDER 
equations for the Fourier coefficients fj, lrn '\p) and 
the LAES amplitude 

1. Results for s-wave scattering (I — 0) 



The low-frequency iterative solution of the integro- 
differential equation (|23|) . presented in Section UTTl cor- 
responds to the iterative account of the integral terms 
M s>s ' in Eq. (|A5|) for solving the system (IA1|) . In the 
lowest order in M SjS >, the solution of Eq. (IA1|) is: 



Equation can be converted into a system of in- 
homogeneous linear algebraic equations for the Fourier 
coefficients /fe(p) of the function / p (r) = J2k fk(p)e~ lkT '■ 

VM s y(e + (5Sw)/ 2sH j(p) = kc 2s+ «5(p), (Al) 



where the symbol d is equal to (1) for an even (odd) 
k. The inhomogeneous term in the system (|A1[) is ex- 
pressed in terms of Fourier coefficients of the wave func- 
tion Xp (r = 0,r) [cf. Eq. JS)]: 

where J n [z, x) is a generalized Bessel function: 



J n (z,x) = J2 e l(n+2p)ars(z) J„ +2 p(M)W- 

p— — oo 

Therefore, the system (|A1|) is equivalent to two separate 
(uncoupled) systems for even and odd Fourier coefficients 
of the QES wave function $ p (r,i) at r — > 0. 

The matrix elements A4 S)S <(e) in Eq. (|A1[) have the 
following form: 



M a , s >(e) = 4 1 (p 2s )^ :S ' - M >ia #(e), 
•A(P2«) = 



M s>8 /(e) = i 



-a 1 + r fc| s /2 - ifc 2s ' 
00 dr 



(A3) 

fc 2s = V -f , (A4) 



x[e- A Wj s _,(fe(T))-4. 



e le »+ 



(A5) 



A(t) = — — (t — snr — 



4 . 2 r 



27' 



/ \ w p / • 4 • 2 T 
z t = r- sm r sin — 

where J n (x) is a Bessel function, and the following no- 
tations are used in Eqs. (|A3|) - (|A5|) : e„ = e + nfiw = 



i? + M p + nfiw, p„ = y/2me n . Note that only diagonal 
matrix elements M s , s ' contain the information on atomic 
dynamics [i.e., the field- free elastic scattering amplitude 
■A(p2s) for a "momentum" p 2s , which is imaginary for 
closed channels, with e 2s < 0], while the non-diagonal 
elements (s 7^ s') depend only on the incident electron 
energy E and the laser parameters. 

In terms of the coefficients /fe(p), the LAES amplitude 
(|26|) can be represented in an alternative form 



A n (p,p n ) = K 1 fk(p)c* k -n(Pn)- 



(A6) 



fk(p) ~ nA(p k )[ck - S,A(pk+2 S >)M 0tS t(ek)ck+2s>] , 

The first term in the approximation (|A7[) corresponds to 
the zero-order approximation (137p for the function / p (r) , 
while the second term describes the rescattering correc- 
tion (|47[) . However, we emphasize that the approxima- 
tion (|A7[) is more accurate than the low-frequency ex- 
pansion (J2U) because the LAES amplitude (|A6I) [as well 
as the sum over s' in Eq. (|AT|> ] involves a summation 
over all intermediate channels, including closed channels. 
Nevertheless, using the approximation (|A7I) we are not 
able to provide a closed-form analytic expression for the 
LAES amplitude. Finally, we note that all non-diagonal 
matrix elements M s , s > (with s 7^ s') are equal to zero for 
a circularly polarized (t — 0) field F(i). In this case the 
sum over s' in Eq. (|A7|) contains only the single term 
with s' = 0. 



2. Results for p-wave scattering (I — 1) 

For 1 = 1, matching the QES wave function (|13|) [with 

$p C \r,t) given by Eq. (fl8|) ] to the small- r boundary 
condition (JTTJ) results in the system of three (for fi = 0, 
±1) coupled integro-differential equations for functions 

jp^V) = E fe /i lM) (p)e" ifcT [cf. Eq. ([231) for the case 
1 = 0]. This system can be converted into the following 



three matrix equations for the Fourier coefficients / ( 
/^ } (P): 



O) 



E 



M s (0 s l(e 5 )/ ; 



(o) _ Ji M 



--2S+5' 



(A8) 



(es) M 





f?f 5 ),(A9) 

■ C 2s+(5 / 



where e$ = e + SHuj and S is equal to (1) for an even 
(odd) k, similarly to the result for s-wave scattering in 

Eq. (|A1|) . The coefficients on the right-hand side of 
Eqs. (|A8|) . (|A9[) can be expressed in terms of the coeffi- 
cients Cfe(p), given by Eq. (|A2|) : 



Jp) 



(P) 



1 



47ry i ;(p)c fc (p) + W3(i+^) 



1 



c fe -i(p)- 1- 



1 + f 



4?kj 
Cfc+i(p) 



k= — 00 



where the spherical harmonic Yi M (p) is defined as in 
Ref. 0. 

The matrix elements M s ( ° } ,(e), M s ( ^,(e) and M s (M j(e) 
(^i = ±1) in Eqs. (|A8|) and (|A9I) have the following form 
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(cf. Ref. Hg): 

Mi5(e)=(-i 

\ Cil 



-C 



dr 

75/2' 



Je s+B ,T/(huj) 



x[e- lX( - T h s - s ,(£z(r))-5 s 



(A10) 



<2 00 +c 



e i<i s + B ,T/(hu))-i\(T) 



x{[ipi(r) + ntz(T)]J s - 3 ,(ez(T)) 
-£p 2 (r)J' s _ s ,(Hr))}, 



(All) 



M { S %)=C 



dr 

7m ( 



i ie s+s> T / (few)— iA(r) 



:{ -^Pi(r)J,_^(fo(T))+ft,(T) [j^ s ,(fe(r)) 



&(t) 



J s — s' 



(£z(r))]} 



(A12) 



where J^(.z) is the derivative of the Bessel function and 
the following notations are used: 



3i s - s ' +1 /mw\3/2 



ft / 



C = 



Pi (t) = — - I —~ sin sin r + cos - 



Mp / 4 . 2 T 2 



Mp / 4 . 2 T 2 



ftw V T 



^( r ) = TT. Z2 sin - - - sinr + 1 



Once the Fourier coefficients fj^\p) are known, the 
exact TDER result for the p-wave LAES amplitude is 
given by: 



-1 fc= 



(A13) 



Appendix B: The uniform asymptotic 
approximation of the integral (|53|) 

In this Appendix, we describe the approach for the 
uniform asymptotic expansion of the integral (|53[) . We 
note first that after replacing the integration variable r 
in Eq. (|53[) by x = t — n/2 — ipt, the amplitude An^ is 
expressed in terms of the integral /„ (p) : 

4°) = jV"*V n (p), 

W) = 7^ f(x)e i ^dx, 



2tt 



(Bl) 



where / (x) = A(x + 7r/2 + (p t ) is a periodic function 
of x and ip(p,x) = nx — pshix. Assuming p > 1 and 
p > |n|, the main contribution to the integral is 
given by the neighborhoods of the saddle points x = x± , 
satisfying the equation dip(x)/dx = 0: 

Ti 

x± = ±a, cosa=— , < a < 7T. (B2) 
P 



Since the points x± tend toward each other and coa- 
lesce at a = 0, following the general idea of the uni- 
form approximations of integrals j3.6j , we rewrite the pre- 
exponential function f(x), explicitly extracting the term, 
which approximates the f{x) in the neighborhood of the 
two coalescing saddle points. Taking into account the 
periodicity of f[x), we rewrite it in the following form: 

f(x) — ao + a\ sin a; + (cos x — cos a)g(x), (B3) 

where ao and a\ are easily determined to be 

f(x+)+f{xJ) f(x+) - f{x J) 
a °- 2 ' ai ~ 2^ ' 

and where g(x) is an analytic, smooth, periodic function 
of x. After substituting Eq. (|B"3)l into Eq. ([Bl]) . the inte- 
gration of the first two terms of the expression (|F33|) can 
be performed analytically. The result for I n (p) is: 

I n {p) = a J„(p) + iaxJ' n {p) + I n (p), (B4) 

where J n (p) and J' n {p) are the Bessel function and its 
derivative, while I n (p) is the remainder integral: 

I„(p) = —[ (cosx-cosa)g(x)e iip ^ x) dx. (B5) 
2tt J-* 

Integrating I n (p) by parts, we obtain 

f n{p) = J- T ^El e Mp^) dx . (B6) 
Amp J_ n dx 

Comparing Eq. (IB6|) with Eq. ()B1|) . one sees that the 
remainder term I n (p) has the same form as the original 
integral (jBll) . but contains a small parameter p~ l . Rep- 
resenting the function dg{x)/dx in Eq. (IB6|) by the form 
(|B3|) and applying the same integration procedure as for 
I n (p), we find the asymptotic expansion of the integral 
I n (p) for the large parameter p. 

For the case of a Kroll- Watson-like approximation, we 
neglect the remainder term I n (p) in Eq. (|B4|) . which 
gives immediately the result (|55[) for the scattering am- 
plitude aI?. 

Also, we recall here another asymptotic approximation 
of the integral (|B1[) , which was suggested in Ref. [34j |. 
where the integration interval in Eq. (|B1[) was divided 
into two parts (— tt > x > and > x > tt) followed by 
taking into account the saddle points x± independently 
(as non-coalescing saddle points). The result is that the 
integral I n {p) can be expressed in terms of the Anger 
function, J n {p), (which coincides with the Bessel function 
for integer n) and the Weber function, E n (p) [371 ]: 



I n (p) = a + J n (p) + ia-E n (p), 
/(*+) ± /(*_) 



(B7) 
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